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Abstract 

We examined the characteristics of a metal-insulator transition on a two-dimensional 
non-bipartite lattice when both disorders and electron interactions are present. Us- 
ing a real-space renormalization group method and finite-size scaling analysis we 
have found that after the disordered system transits from the insulating state to 
the metal one, it can reenter an insulating state if there is a further increase in the 
electron-electron interactions. 
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Ever since the discovery of the metallic properties on effectively 2D electron 
systems in metal-oxide-semiconductor field-effect transistor (MOSEF's) [f], 
there have appeared a spate of discussions and investigations over its underly- 
ing physical mechanism [2] . In experiments, as the carrier density n is increased 
above a critical density ric, the conductivity adc decreases upon lowering the 
temperature T (typical of an insulator) and as n < ric, CFdc increases upon 
lowering T (typical of a conductor). Because the data can be scaled onto two 
curves, one for the metal and the other for the insulator, this phenomenon is 
seen as evidence for the occurrence of a quantum-phase transition, namely the 
metal-insulator transition (MIT), with carrier density as the tuning parame- 
ter. But there are also other viewpoints, supported by some experiments [3,4], 
maintaining that there is no 2D MIT at zero temperature and the metallic 
behavior is attributed to the conventional, though nontrivial, electron trans- 
port. For example, two of these models are based upon electron scattering by 
impurities [5] and the effect of temperature dependent screening [6,7]. Thus 
there is, as yet, no general consensus about the origin of the metallic behavior 
and it remains a controversial topic. 
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According to the conventional one-parameter scaling theory of noninteracting 
electrons [8], any amount of disorders will localize the 2-D electronic state 
and make it insulating. Moreover, the metallic state has been seen in highly 
mobility structures at low carrier concentrations. Thus it is expected that the 
combining effect of both electron-electron interactions and the disorders are 
very likely to play an important role. In fact, the role of electron-electron inter- 
actions in disordered systems was recognized long ago by Finkclstcin [9] and 
Castellanic [10], who treated the disorder in lowest order, where all interac- 
tions contributed t o the leading logarithmic behavior were summed. Recently, 
Si and Varma [11] calculated a correction to the compressibility of a disordered 
system by considering the ring diagrams. Because of the perturbation nature 
of these methods, they can not deal with the regime with the disorder and 
the interactions having comparable magnitudes. Hence most recent work has 
used exact diagonahzation [12,13] and Monte Carlo methods [14] in the inves- 
tigations. As we know, both methods suffer from intensive calculations and it 
is very difficult to apply them to large-size systems. To resolve this difficulty, 
we resort to the real-space block renormalization group (BRG) method[19]. 
Although this method has uncontrollable approximations, it can give us many 
qualitative and insightful results and show us the direction for further more 
accurate work. The model we use is the Anderson-Hubbard model, which is 
the most simple one to include the essential ingredients for our purpose, i.e. 
the disorder and the electron-electron interactions. We will use this model to 
study the charge gap and the electron localization and delocalization on a 
triangular half-filled lattice as shown in Fig.l. Normally there are two kinds of 
disorders. One is the site disorder related to the site fluctuations and the other 
is the bond disorder related to the hopping terms. Here we will only consider 
the former one. The Hamiltonian for Anderson-Hubbard model in our case 
can be written as, 
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where t is the nearest-neighbor hopping (exchange coupling) term, U is the 
local repulsive interaction and /i is the chemical potential, c^(cjcr) creates 
(annihilates) an electron with spin a in the valence orbital of the dot located 
at site i; the corresponding number operator is riia = cf^Ci^j. < • • • > on the 
first sum in Eq. (1) indicates that summation is restricted to nearest-neighbor 
dots. H.c. denotes the Hermitian conjugate. Note that this model Hamiltonian 
allows only one orbital per dot. That orbital can be empty or accommodate 
one or two electrons. U is the repulsion of two electrons (of opposite spins) 
placed in the same dot. Si measures the fluctuation of the site energies. It is 
assumed that Si follows a Gaussian distribution with the width to be W, i.e.. 
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in which the bar over e means its average value. Initially, we use £ = 0. 

The essence of the BRG method is to map the above many-particle Hamil- 
tonian on a lattice to a new one with fewer degrees of freedom and with the 
same low-lying energy levels [18]. Then the mapping is repeated leading to a 
final Hamiltonian of a seven-site hexagonal array for which we obtain an exact 
numerical solution. 

When there is no disorder, namely Si — 0, the procedure can be summa- 
rized into three steps: First divide the A^-site lattice into appropriate n^-site 
blocks labeled by p (p=l,2,..., N/ug) and separate the Hamiltonian H into a 
intrablock part Hb and an interblock Hjb, 

H^Hb + Hib^ Y.Hp + E (3) 

p {p,p') 

where Hp is the Hamiltonian (1) for a given block and the interblock p,p' 
coupling is defined in Eq. (4) below. 

The second step is to solve Hp exactly for the eigenvalues Ep^ and cigcnfunc- 
tions $pj(i = 1,2, ...,4"^). Then the eigcnfunctions of Hb arc constructed by 
direct multiplication of . The last step is to treat each block as one site 
on a new lattice and the correlations between blocks as hopping interactions. 

The original Hilbcrt space has four states per site. If we are only concerned 
with lower lying states of the system as when studying the metal-insulator- 
transition [17], it is not necessary to keep all the states for a block. 

To make the new Hamiltonian tractable, the reduction in size should not 
be accompanied by a proliferation of new couplings. Then one can use an 
iteration procedure to solve the model. To achieve this, it is necessary to keep 
only 4 states in step 2. Their energies are Ei (z=l,2,3,4). In order to avoid 
proliferation of additional couplings in the new Hamiltonian, the four states 
kept from the block cannot be arbitrarily chosen. Some definite conditions as 
discussed in Ref. [15] must be satisfied. For example, the states must belong 
to the same irreducible representation of Cqu symmetry group of the lattice. 
In particular, in order to copy the intrasite structure of the old Hamiltonian, 
Es = i?4 is a necessary condition. Furthermore, particlc-liolc symmetry of a 
half-filled lattice requires that Ei = E2. Further restrictions follow from the 
need to make extra couplings vanish. Operators in the truncated basis are 
denoted by a prime so that the interblock coupling of Eq. (2) is 

^pp' = ^A2^E4Xv> (4) 
where v represents the number of couplings between neighboring blocks. The 
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coupling strength for the border sites of a block by A and the renormalization 
group equation for the coupling strength is 

t' = v\H, (5) 

The other renormalization relation is 

U' = 2{Ei-E2). (6) 

Once we introduce the disorders, because there is no exact particle-hole sym- 
metry any more, the parameters t, U and have to be renormalized on aver- 
age. In details, let us use a and /? to be the block indices. Then for one block 
a, we can have 



U'^^E^ + -2E^, (7) 
e'^^E^-E^. (8) 

After the renormalization, the new energies do not obey a Gaussian dis- 
tribution. In order to iterate the RG, as in Ref.[19], we adopt the following 
procedures: 

(1) is forced back into a Gaussian distribution with the new width 



W' = (£")2 -(i^f. (9) 



The new Gaussian is not centered at zero since there will be a constant shift 
due to the electron interactions. But we can still take it to be zero by formally 
introducing the chemical potential. 

(2) 17°' is forced back to be constant. 

U' = W. (10) 

(3) To get the renormalized hopping term, wc will have to consider all the 
possible non-zero average values of the coupling between the block states. For 
two neighboring blocks, there are 4 possibilities. 

tf=t<ro^i E [cL.,,>c/.,.]if o">, (11) 

<ai,Pj>,a 

tf=t<Tro^l E f r>, (12) 
tf=t<ri^\ E [d.%.]in''o">, (13) 

<ai,l3j>,a 



4 



tf=t<Tn^i E Kiacpj.\\nn''> . (14) 

<ai,Pj>,a 

We also force the distribution of into a Gaussian with mean 

t' = If, (15) 

and width 

t', = J^-(f)\ (16) 

In principle, to use more blocks for averaging is always desirable. In our work, 
we use 7 unit blocks, which themselves form a hexagonal block. 

For clean systems, we have shown in our former work [15] the existence of 
Mott MIT at f/ = 12.5 by investigating the charge gap Ap, which can be 
expresses as the limiting value of the renormalized U, 

Ag^lirrir,^^ (17) 



For disordered systems, as we know from the orthodox scaling theory [8], 
the nonintcr acting electrons will be localized, which leads to another kind 
of insulator, namely the Anderson insulator. The localization feature can be 
exposed by the inverse participation rate (IPR) [16] ^ defined as 

^(4) ^ ^ _^ (18) 

E \ni\ 



where Ng is the total number of the electrons, which is equal to the site number 
Ng and ip the system wave function. For totally localized electronic states in 
the half-filled system, half of the sites have (if) |nj| ip) = 2. On the other hand, 
when the electron is totally delocalized, the average number of electrons per 
site should be {^p \ni\ ip) = 1. Hence the IPR will satisfy 0.125 < < 1 
between the above two limiting cases. 

In Fig. 2, we demonstrated how the disorders infiucncc the Mott MIT. It is 
interesting to note that the disorder shows the effect of destablizing the Mott 
insulating state, namely, as we increase the disorder, part of the insulating 
region appears "melted" into a "metallic" state, in which the charge gap is 
nearly zero. At first sight, this feature quite contradicts the generally accepted 
viewpoints that the disorder should stabilize the insulating state. But actually 
they are talking about different interaction regimes with two different MIT 
mechanism. Our result is obtained in a strong coupling regime with the charge 
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gap to measure MIT (Mott MIT). As will be shown below, for weak electron- 
electron interactions, the disorders really stabilize the insulating state with 
^("^^ to signature MIT (Anderson MIT). 

Fig. 3 displays the evolving of electron localization and delocalization as we 

tune both the disorder and electron-electron interaction in the whole coupling 
regime. It is easy to see that all the curves follow almost the same variation 
patterns as we increase U for fixed disorder W. Generally speaking, there are 
three stages. 

1) When U is very small, the electron-electron interactions will help with 
the electron delocalization. This is physically understandable since the on- 
site repulsion will tend to rearrange the charges and make the charge-density 
more homogeneous. This result is consistent with the findings by Ma [19] and 
Caldara [20], but is different from the results from Berkovits [21]and Benenti 
[22]. The reason for this difference might be because Berkovits and Benenti 
use spinless electrons and only consider the neighbor-site repulsion, which does 
not have the homogenizing effect like the on-site interaction. 

2) As we increase U further to about U ~ W. the electron become maximally 
delocalized and becomes localized again when U > W. The delocalization 
peak in the diagram is very apparent when W is, again when U is, for ex- 
ample, bigger than 10. If we interpret the shallow region around the peak 
to be in a metallic state, the system will then experience here a three-phase 
transition: insulator-metal-insulator, which is quite consistent with the exper- 
imental findings [24,25]. Another interesting fact we might note is that, as the 
disorders increases, the metallic region moves further to the right of the axis, 
which implies that the initial insulating state becomes more stabilized. This 
is what we have mentioned before. 

3) After the two phase transitions, the system continues to exhibit Anderson 
insulating behavior until finally it evolves into a Mott insulating state with 
very strong electron couplings. In this regime, every site has one electron on 
average because of the strong on-site electron interactions. 

From the above discussion, we should be able to summarize four different elec- 
tronic features as U increases from zero to infinity, namely, Anderson Insulator- 
Metal-Anderson Insulator-Mott Insulator. This full scenario from weak to 
strong coupling is quite consistent with the present experiments [24,25] . 

One of the advantages of the BRG method is that it can be adapted easily to 
carry out finite-size scaling analysis of the system [23] . Thus more fundamental 
physics can be explored and displayed. 

Fig. 4 presents our calculation results with respect to the charge gap for each 
fixed disorder. It is very interesting to note that as we increase the disorders. 
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the finite-size scaling analysis provide us one critical point while W is small, 
which clearly represents quantum phase transition. Then the system enters 
some disorder/interaction range where no clear critical transition point can 
be identified. We call this region a "mixed" region since the two pure phases 
that is necessary to characterize the phase transition are difficult to define, 
this is demonstrated by the absence of a clear crossing point from the finite- 
size scaling analysis. As the disorders become stronger, two critical points are 
evolved out of the mixed region. From here, we can easily identify the two 
phase transitions as discussed above, i.e. insulator-metal and metal-insulator 
transitions. But the last phase transition, Anderson insulator- Mott insulator 
is still absent in the finite-size scaling diagrams up to the magnitude of the dis- 
orders we have considered. It should be expected that one more critical point 
might appear if stronger disorders are considered. More extensive calculations 
are needed before a final conclusion can be drawn. 

In summary, by using the BRG technique, we have investigated in detail the 
delicate interplay between disorders and electron-electron interactions on a 
half-filled triangular quantum dot lattice. The insulator-metal-insulator tran- 
sition has been well demonstrated by the finite-size scaling analysis. 
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Figure 1 




Fig. 1. Schematic diagram of the triangular lattice with hexagonal blocks. Only 
two neighboring blocks p and p' are drawn here. The dotted lines represent the 
interblock interactions and the solid lines the intrablock couplings. 



Figure 2 




Fig. 2. Influence of the disorder (measure by W/t) over the charge gap A^/t depen- 
dence upon the electron-electron interactions (U/t). 
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Figure 3 
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Fig. 3. Variation of the inverse participation rate (IPR, S,^'^^) against the elec- 
tron-electron interactions for different disorders with (a) W/t=l (solid line), 
3(dashed line) and (b) W/t=5(solid line), 10(dashed line). 
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Figure 4 




Fig. 4. Finite-size scaling analysis over the charge gap A^/t against electron-electron 
interactions for fixed disorders, (a) W/t=l. (b) W/t=5. (c) W/t=10. 
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